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Abstract 

We analyze a random resistor-inductor-capacitor (RLC) lattice model of 2 di- 
mensional metal-insulator composites. The results are compared with Bruggeman's 
and Landauer's Effective Medium Approximations predictions where a discrepancy 
was observed for some frequency zones. Such a discrepancy is mainly caused by the 
strong conductivity fluctuations. Indeed, a two-branches distribution is observed 
for low frequencies. We show also that at p c the so-called Drude peak vanish by 
increasing the system size and increase for vanishing losses. 

Keywords: conductivity, composites, resistor networks, percolation. 

1 Introduction 

The optical properties of composite materials are very different from those of their 
constituents particularly in metal insulator composites [U EJ. Theses differences 
are mainly observed at percolation threshold p c at which one of the two compo- 
nents of the composite first forms a connected path extending through the sys- 
tem Bruggeman's and Landauer's Effective Medium Aproximation (BEMA (4], 
LEMA 0) provides the most useful way to describe these properties. 

Zeng et al. (H), following the work of Koss and Stroud |2j, analyzed the optical 
properties of a two-dimensional normal-metal-insulator composite. 

They found a surface plasmon resonance peak of the real part of the conductivity 
for various concentrations above and below p c in agreement with LEMA. For dc 
case, Watson and Leath |H] gave equation to correct the slope of the curve at the 
critical concentration. Koss and Stroud generalized this equation to their ac prob- 
lem. The so-called surface plasmon modes are known to be an important absorption 
mechanism in random metal-insulator composites jH]. Drude peak, corresponding 
to nonzero dc conductivity at and above p c was also found. This peak is due at p c to 
the finite size effects and is expected by Koss and Stroud to vanish in the limit of in- 
finite systems. For system size 100 x 100 Zeng chose a characteristic relaxation time 
value of 10 in RLC systems corresponding to a large loss in the metallic component 
R = 10 _1 (r = L/R). For vanishing losses the local field and the ac conductivity 
fluctuations are strong mainly at percolation threshold p c [THIEI]. Thus calculations 
for only 5 samples, as done by Zeng, for vanishing losses, will not be statistically 
sufficient. 
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Zeng et al. used Frank and Lobb algorithm [12], for frequencies ranged from 1CT 2 
to 1 (expressed in units of the plasmon frequency u p ) . This algorithm is very fast 
but based on many division operations, which causes a numerical instabilities both 
for smaller frequencies and losses. In order to compute the effective conductivity 
for the corresponding 2d resistor network for smaller frequencies and losses, we use 
a method of exactly solving Kirchoff's equations, extensively described in Refs [1 3j . 
By using this method the current is preserved in all nodes of the square lattice and 
the total entering and outgoing currents of the sample are accurately identical. 

In this paper we present a numerical study on 2d resistor network model for a large 
frequency range including the plasmon frequency u p around the bond percolation 
threshold p c for smaller loss and larger size. 

Confirming Koss and Stroud prediction we show that the Drude peak decreases 
by increasing the system size. We show also that for vanishing losses Drude peak 
increase and, at percolation threshold p c , for low frequencies the distribution of the 
real part of the effective conductivity is splitted in two branches, the first one is 
due to the non percolate samples and the second branch to the percolate ones. The 
results allow us to test the validity of the above cited BEMA and LEMA. 

2 Model description 

It is well known that in the optical frequency region a semi-continuous metallic film 
can be modelled as a 2d RL — C lattice [2j. Capacitors C stand for the dielectric 
grains that have dielectric constant e^, while the inductance L represents the metallic 
grains characterized for Drude metal by the following dielectric function [Tl] . 



where e& is the contribution to e m due to the inter-band transitions, u T = 1/r is 
the relaxation rate and uj v is the plasmon frequency. While for dielectric grains, the 
dielectric function is, 



The metal conductivity a m = —ie m Loj^K is characterized by positive imaginary 
part in the optical frequency region (inductive behavior) [2]. We can then model the 
metallic grains as inductances L with loss R (a m = (R+iLu)" 1 ) while the dielectric 
grains can be represented by capacitances C (cr^ = iCuj). We can take without loss 
of generality, L = C = uj p = 1, (uj p corresponds to the case where | a m |~| ad \ 
for vanishing losses, i.e., Cuj p = 1/Lu p ). The inductance L and the capacitance C 
can be assumed to be constant in the range between 1/r and u> p , we express in this 
paper the frequency in units of the plasmon frequency (Co = uj/uo p < 1). 

We model the metal dielectric composite films as a square lattice composed of 
metallic bonds with conductivity : 



e m (uj) = e b - (u} p /u) 2 /[l + iu T /u] 



(1) 



e d = cte 



(2) 



(R + iu) 



-i 



(3) 



and concentration p, and dielectric bonds with conductivity: 



(4) 



and concentration 1 — p. 

In order to compare with Zeng's results, we model the metal by a capacitance 
parallel to RL branch, this is called RL\\C model to be distinguished from the RL 
one. 

3 Results and discussion 

We consider a lattice of size nxn (n = 100, 200 and 600). If we choose R = 10~ 3 , the 
smaller system size is larger than the correlation length (following Zekri et al. [llj, 
the correlation length behaves as R~ l l 2 ). We choose then losses R not smaller than 
the above value so that we ensure a correlation length always smaller than the 
system size. This avoids the size dependent effects of the conductivity. The effective 
conductivity here is averaged from 550 samples and compared to Zeng's simulations 
for only 5 samples. 

Let us first reproduce Zeng's results on the effective conductivity by exactly solving 
Kirchoff laws [THJ for an RL\\C model. We compare them to RL model at p c with 
R = 10" 1 . In Figurela, the two models seem to produce identical results except 
near uj = 1, where the real part of the conductivity seems to vanish for RL\\C model 
while it saturates for RL model, the divergence between the two models start from 
the frequency uj = 0.35. Large fluctuations of this conductivity are observed below 
uj = 10~ 2 . We show also that the two models are in a good agreement with BEMA 
for all the frequency range investigated and are different from LEMA. Comparison 
between the two models at p c with R = 10~ 3 is done for sizes 100 and 200. They 
produce also identical results but the fluctuations extends to the whole frequency 
range investigated. Then we have done calculations for large samples using RL 
model. In Figurelb comparison between RL model for system sizes 100, 200 and 
600 and both BEMA and LEMA is done at p c with R = 10~ 3 . The size effect on 
Drude peak is also investigated. The results obtained and BEMA are alike only 
for frequencies ranged from 10~ 2 to 1, and they follow qualitatively LEMA for low 
frequencies, where Drude peak increase by decreasing the loss. We see also in this 
figure that Drude peak decreases by increasing the system size and should vanish 
in the limit of infinite systems as predicted by Koss and Stroud 0. The data in 
the described curves are directly averaged (i. e. Arithmetic Average (AA) of the 
550 samples for sizes 100 and 200 and 30 samples for the size 600). In order to 
check the validity of this average, the same results for the size 200 are plotted using 
Logarithmic Average (LA). The LA curve and the AA one are alike for frequencies 
ranged from 10~ 2 to 1, and they diverge for low frequencies. 

It is then important to investigate the conductivity distributions for frequencies 
in the range examined here. For u> = 1 the distribution of the real part is Gaussian 
centered at 1. For uj = lO^ 1 the distribution becomes broader and at uj — 10~ 2 
it becomes poissonian (the distributions for uj = 1, 10 _1 and 10~ 2 are not plotted 
here). In Figures2a we show for uj = 10~ 3 that the distribution is splitted in 
two branches, the first one has the form of a peak centered at the vicinity of zero 
and the second one has a low height and is broad, its width tends towards large 
positive values. The two branches distribution observed for frequencies less than 
10~ 2 are the reason of the large divergence between AA and LA. In Figures2b and 



2c we compare for the frequency Q = 1CT 3 the distributions of the real part for 
concentrations p of the metallic grains (RL bound) below and above p c . Above p c 
means p—p c = and below p c , p—p c = —1/N 1 ^ where the connectivity length 

exponent v [2] is 4/3 for d = 2 and N is the system size. The peak centered in the 
vicinity of zero (first branch) increases below p c and tends to disappear above p c , 
whereas the width of the second branch tends to disappear below p c , and increase 
above p c and becomes dominant. Therefore the two branches appear only at p c . 
The first branch is related to the non percolate samples and the second branch is 
related to the percolate ones. Similar behavior of this conductivity distribution has 
been observed recently [15j for the distribution of the critical links (singly connected 
links pi El). The disappearance of one branch corresponds to that observed for 
critical links below p c (also observed when the system size increase (Hj). Such 
branch corresponds to configurations for which the infinite cluster has a very small 
number of critical links. 

4 Conclusion 

In this paper we have compared (RLC) lattice model with BEMA and LEMA at and 
near the percolation threshold p c for large and small loss in the metallic component. 
For large loss the results and BEMA are alike for all the frequency range investigated, 
but for small loss they are alike only for frequencies ranged from 10~ 2 to 1, and they 
follow qualitatively Landauer's approximation for lower frequencies. We have also 
shown that for small loss the arithmetic and the logarithmic average of the real 
part of the effective conductivity diverge for frequencies less than 10~ 2 due to two 
branches distribution of the real part of the effective conductivity. The first one is 
due to non percolate samples and the second branch is due to percolate ones. We 
showed also in this paper that Drude peak increases by decreasing the metallic loss 
and vanish with the size as predicted by Stroud fZj at p c . 
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Figure 1: a) R = 1CT 1 , Arithmetic Average (AA) over 5 samples of size 100. b) 
R = 10~ 3 , RL model, AA over 550 samples for the sizes 100 and 200 and 30 samples 
for the size 600. Logarithmic Average(LA) for the size 200. The real part of the 
effective conductivity is expressed in arbitrary units and the frequency is expressed 
in units of the plasmon frequency. 
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Figure 2: Distribution of the real part of the effective conductivity expressed in 
arbitrary units, using RL model and AA over 550 samples of size 200 at frequency 
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